DOE/ER/03077-228 


■•-I  en 
13  >J  C 

>H  O 

2        O 


Courant  Mathematics  and 
Computing  Laboratory 


U.S.  Department  of  Energy 


Construction  of  Solutions  for 

Two  Dimensional  Riemann  Problems 

W.  B.  Lindquist 


CM 

w 

• 

c 

o 

o 

•H 

4-> 

00 

a  rH 

Ol 

.-H  (0 

CM 
1 

s! 

•H 

r~ 

r~ 

y-i  w 

o 

o  c 

ro 

0) 

o 

C  § 

\ 

4-)  O  •'-' 

« 

tn  -H  ^3 

w 

•:!^  ^ 

Research  and  Development  Report 

Supported  by  the  Applied  Mathematical  Sciences 

subprogram  of  the  Office  of  Energy  Research, 

U.S.  Department  of  Energy  under  Contract  DE-AC02-76ER03077 

°  -S  m  1-.        Mathematics  and  Computers 


September  1984 


New  York  University 


UNCLASSIFIED 


DOE/ER/03077-228 

UC-32 

Matnematics  and  Computers 


Courant  Mathematics  and  Computing  Laboratory 
New  York  University 


CONSTRUCTION  OF  SOLUTIONS  FOR  TWO  DIMENSIONAL  RIEMANN  PROBLEMS 

W.  B.  Lindquist 
September  1984 


Supported  by  the  Applied  Mathematical  Subprogram 

of  the  Office  of  Energy  Research, 

U.  S.  Department  of  Energy,  Contract  DE-AC02-76ER03077 , 


UNCLASSIFIED 


-11- 


DISCLAIMER 


This  report  was  prepared  as  an"  account  of  work  sponsored 
by  an  ap;ency  of  the  United  States  Government.   Neither 
the  United  States  Government  nor  any  aerency  thereof,  nor 
any  of  their  employees,  makes  any  warranty,  express  or 
Implied,  or  assumes  any  lepal  liability  or  responsibility 
for  the  accuracy,  completeness,  or  usefulness  of  any 
Information,  apparatus,  product,  or  process  disclosed,  or 
represents  that  Its  use  would  not  lnfrinp;e  privately 
owned  rights.   Reference  herein  to  any  specific  commercial 
product,  process,  or  service  by  trade  name,  trademark, 
manufacturer,  or  otherwise,  does  not  necessarily  constitute 
or  Imply  Its  endorsement,  recommendation,  or  favoring  by 
the  United  States  Government  or  any  agency  thereof.   The 
views  and  opinions  of  authors  expressed  herein  do  not 
necessarily  state  or  reflect  those  of  the  United  States 
Government  or  any  agency  thereof. 


Printed  in  U.S.A. 

Available  from 

National  Technical  Information  Service 
U.S.  Department  of  Commerce 
5285  Port  Royal  Road 
Springfield,  VA  22l6l 


2, 


-Ill- 


Table  of  Contents 


Abstract  1 

1.  Introduction  2 

2.  The  2-Dimensional  Riemann  Problem  3 

2.1.  The  Two  Dimensional  Nonlinear  Waves  for 

Piecewise  Smooth  Solutions  5 

3.  Two  Dimension  Merhod  of  Solution  Construction  - 
Examples  11 

3.1.  A  Convex  f  Example  12 

3.2.  A  Single  Inflection  Point  Example  13 

4.  The  Two  Dimensional  Riemann  Problem  for  Two  Phase  Flow 

in  Porous  Media  15 

5.  Concluding  Remarks  17 
Acknowledgments  17 
References  18 
Figure  Captions  19 
Figures  21 


Construction  of  Solutions  for  Two  Dimensional  Riemann  Problems 

W.  B.  Lindquist " 

Courant  Institute  of  Mathematical  Sciences 

New  York  University 

New  York,  N.Y.  10012 


ABSTRACT 


Solutions  to  the  scalar  quasilinear  equation 

forfiiC^-.R-'R  with  initial  data  given  by  a  two  dimensional  Riemann  problem  are 
piecewise  smooth  if  fi"  fi^  f,  and  /  has  at  most  one  inflection  point.  We  show  the 
'pieces'  of  this  solution  can  be  classified  and  are  expressible  in  terms  of  two  dimensional 
non-linear  waves  in  analogy  with  the  non-linear  rarefaction  and  shock  waves  of  the 
Riemann  problem  in  one  spatial  dimension.  The  two  dimensional  waves  can  be  expressed 
in  almost  closed  form.  Explicit  solutions  are  constructablc  from  these  waves.  An 
application  is  illustrated  by  calculation  of  the  interaction  of  water/oil  banks  in  two  phase 
incompressible  flow  in  reservoirs. 


a.  Supponed  by  the  Applied  Mathematical  Sciences  subprogram  of  the  Office  of  Energy  Research, 
U.  S.  Dept.  of  Energy,  Contraa  DE-AC02-76ERO3077. 
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1.  Introduction 

The  Cauchy  problem 

Mr^  ^  J-^t^MLf)!  =  0  ,  (1.1) 

at        ^.1      dx, 

with  initial  data  piecewise  constant  on  a  finite  number  of  wedges  focused  on  a  single  point  in  the 
X,  y  plane  is  deflned  as  a  two  dimensional  Rienumn  problem.  Without  loss  of  generality,  this  point 
can  be  taken  to  be  the  point  x  =  0,  y  =  0.  In  a  previous  paper  [1]  we  have  shown  that,  for 
ff-.R-R'm  C^,  if  /i  ■  /2  ■  /,  with / having  at  most  one  inflection  point,  the  solution  to  the  two 
dimensional  Riemann  problem  is  piecewise  smooth.  The  proof  consists  of  two  parts.  First  it  is 
shown  that  the  two  dimensional  Riemann  problem  for  /^  ■  /2  ■■  /  is  equivalent  to  a  generalization 
of  the  Riemann  problem  in  one  dimension.  It  is  then  proven  that  the  generalized  one  dimensional 
Riemann  problem  is  piecewise  smooth  if  /  has  at  most  one  inflection  point.  The  proof  of  the 
smoothness  of  the  solutions  to  the  two  dimensional  Riemann  problem  does  not  provide  a 
convenient  method  of  constructing  its  solutions.  In  this  paper  we  formulate  a  construction  that  is 
purely  two  dimensional  and  has  the  advantage  (at  least  in  specific  cases)  of  generalizing  to  the 
problem  /j  =^  /j.  This  construction  consists  of  identifying  the  general  two-dimensional  non-linear 
waves,  analogous  to  the  rarefaction  and  shock  waves  of  the  one  dimensional  Riranann  problran, 
and  then  using  the  method  initiated  by  Guckenheimer  [2]  and  Wagner  [3]  to  piece  these  waves 
together  into  an  entropy  obeying  solution. 

Existence  and  uniqueness  of  solutions  to  the  two  dimensional  Riemann  problem  within  the 
dass  of  bounded,  measurable  functions  is  due  to  Kruzkov  [4].  For  piecewise  smooth  solutions, 
Kruzkov's  uniqueness  condition  reduces  to  two  requirranents  on  the  jump  discontinuities  in  the 
solution: 

n[«*-«-,  /(O -/(«-)]  =  0,  (1.2) 

and 

n[k-u-,Ak)-Au-)]^0,  (1.3) 

where/  stands  for  the  vector  [/i,/2  ]•  In  (1.2),  the  normal  n  to  the  discontinuity  is  oriented  such 
that  u~  s  u*  and  Jfc  is  any  constant  such  that  u~  s  it  «  «*.   We  refer  to  [1]  for  definition  of  the 


terms  used  here  and  throughout  this  article. 

Equation  (1.2)  is  denoted  the  jump  condition  for  the  discontinuity  and  (1.3)  the  entropy 
condition;  in  /?'  x  R'^  (1.2)  is  the  familiar  Rankine-F&igoniot  condition  and  (1.3)  is  equivalent  to 
the  entropy  condition  of  Oleinik  [5].  The  local  nature  of  conditions  (1.2)  and  (1.3)  allow 
construction  of  piecewise  smooth  global  solutions  by  piecing  together  local  solutions  which 
individually  obey  (1.2)  and  (1.3). 

In  §2,  the  two  dimensional  nonlinear  waves  are  described  for  the  class  of  problans  /j  ■  /j. 
In  §3  a  method  for  explicitly  constructing  the  2-dimensional  solutions  in  this  class  from  the 
nonlinear  waves  is  described.  Complete  solutions  for  2-dimensional  Rionann  problems  for 
representative  forms  of  the  function  /  are  given.  Solutions  to  the  problem  of  water/oil  bank 
interactions  in  two  phase,  incompressible,  gravity  free,  flow  in  reservoirs  can  be  obtained  from  the 
study  of  these  2-dimensional  Riemann  problems.  The  results  for  this  physical  problem  are 
presented  in  §4.  Conclusions  and  conjectures  for  the  solutions  to  2-dimensional  Riemann 
problems  for  the  case/  having  more  than  one  inflection  points  or/^  #  /j  are  presented  in  55. 

2.  The  2-Dimeiisioiial  Riemaiio  ProMem 

For  convenience  we  shall  use  the  notation /j  -  /,  /j  —  g,  x^  ■  or,  jcj  "  y. 

A  sufficient  condition  for  a  two  dimensional  Riemann  problem  to  have  a  picewise  smooth 
solution  is  given  by  [1] 

Theorem  2.1:  The  imique  (in  the  sense  of  Kruzkov)  solution  in  the  plane  y  >  x,  to  (1.1) 
with  initial  data  that  is  piecewise  constant  on  a  finite  number  of  wedges  focused  on  a  single  point 
in  the  plane,  with/,  -  /j  -  /,  /  €  C^  :  /?  -  /?, /  having  at  most  one  inflection  point,  is 

a)  piecewise  smooth. 

b)  composed  of  non-linear  waves  and  constant  states  v/idch  are  the  images  under  a  continuous 
map  M  of  the  rarefaction  and  shock  waves  and  constant  states  of  a  generalized  Riemann  problem 
in  one  dimension. 

c)  composed  of  curves  of  irregular  points  corresponding  to  images  under  the  map  M  of  the 
irregular  points  in  the  1 -dimensional  Rionann  problon. 


For  the  case  g  ■  /,  under  the  45°  rotation  2^  =  x  +  y,  2^  =  y  -  x,  (1.1)  becomes 

ilL  +  MA  =  0  .  (2.1) 

From  (2.1)  we  see  that  the  solution  can  be  obtained  along  each  t|  =  constant  plane  independent  of 
other  Ti.  In  particular,  this  implies  the  solutions  obtained  in  the  t)  <  0  half  space  can  be  obtained 
independently  of  the  solutions  in  the  -q  >  0  half;  by  the  symmetry  of  the  problem,  no  new 
solutions  will  be  foimd  in  the  t|  <  0  half  space  that  are  not  found  in  the  upper.  We  therefore 
restrict  our  discussion  to  the  half  space  t]  >  0  (y  >  x).  The  solution  in  the  plane  ti  =  0  is  given 
in  Lemma  2.2  below.  (See  [1]) 

Lemma  2.2: 

a)  K  the  half  line  t]  =  0,  ^  <  0  is  not  a  line  of  discontinuity  of  the  initial  data  (ie:  is  not  a 
wedge  line)  the  waves  incident  upon  the  corresponding  half  plane  f  >  0,  ^  <  0,  -ri  =  0  are 
continuous  across  the  half  plane. 

b)  If  the  half  line  ti  =  0,  ^  <  0  is  a  line  of  discontinuity  (constant  jump)  of  the  initial  data,  it 
remains  a  half  plane  r  >  0,  ^  <  0,  ti  =  0  of  (in  general  variable)  jump  discontinuity  in  the 
solution. 

Similar  statements  hold  for  the  half  line  ti  =  0,  ^  >  0. 

Although  the  map  M  mentioned  in  Theorem  2.1  provides  a  means  of  constructing  the 
solution  to  the  2-dimensional  Riemann  problan,  it  is  not  a  convenient  method  of  doing  so.  In 
particular  it  is  a  method  that  will  not  generalize  to  cases  f  =^  g  where  no  appeal  to  a  one 
dimensional  analysis  can  be  made.  We  therefore  proceed  by  formulating  a  construction  that  is 
purely  two  dimensional  and  has  the  advantage  of  generalizing  (at  least  in  specific  cases)  to  the 
problem  f  =^  g.  We  propose  that  the  correct  method  of  dealing  with  the  general  solution  to  the 
2-dimensional  Riemann  problem  is  to  identify  the  general  two  dimensional  non-linear  waves.  This 
construction  method  (initiated  by  Guckenheimer  [2]  and  Wagner  [3]  )  can  then  be  used  to  piece 
these  waves  together  into  an  entropy  obeying  solution  in  a  manner  that  places  no  reliance  on  one 
dimensional  analyses. 
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2.1.  The  Two  Dimensional  Non-Linear  Waves  for  Piecewise  Smooth  Solutions 

We  proceed  with  the  analysis  and  definition  of  rarefaction  and  shock  waves  valid  for 
piecewise  smooth,  entropy  obeying  solutions  to  the  2-dimensional  Riemann  problem  for  the  case 
/  ■  f .  The  explicit  forms  of  these  waves  and  the  complete  characterization  of  the  irregular  points 
are  given. 

The  2-dimensional  Riemann  problem  is  invariant  under  the  similarity  transformation 
(^  X,  y)  -  (ct,  ex,  cy).  The  solution  is  therefore  constant  along  rays  having  the  space-time  origin 
as  vertex.  Consequently  the  solution  can  be  determined  by  its  restriction  to  any  plane 
t  =  constant(>  0).  The  plane  /  =  1  is  most  convenient  for  this  purpose. 

In  addition  to  this  self-similarity  property,  at  points  at  which  u  is  regular,  the  solution  is 
constant  on  characteristics  (f,  x(t),  y(t))  which  are  the  straight  lines 

—rr  =  ^'(")  •     1-^  =  «'(")  .'0  =  0.  •  (2.2) 

Shock  Waves  and  the  Entropy  Condition 

Let  (t,  X,  y)  be  a  point  of  jump  of  a  unique  piecewise  smooth  solution  to  a  2-dimensional 
Riemann  problem.  Condition  (1.2)  thus  states  that  the  vector 

[l,5/«^«-),v^«")]. 

where 

is  a  tangent  vector  at  the  point  of  jump  t,x,y.   However,  by  the  self -similarity  of  the  Riemann 
solution,  the  vector 

[  1,  x/t,  yit  ] 
must  also  be  tangent  to  the  point  of  discontinuity.  Consequently,  the  vector 

[  0,  x/r  -  5/«^«-),  yit  -  5/u^«-)  ]  (2.3a) 

is  tangent  to  the  point  of  discontinuity.  In  the  plane  f  =  1,  (2.3a)  has  the  simple  form 
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[  0,  X  -  5/u*,«-).  y  -  5,(«+,«-)  ]  .  (2.3b) 

Let  a  >  b.  Consider  all  planes  passing  through  the  points  (1,  Sy(a,b),  Sg(a,b))  and  the 
origin  (0,  0,  0).  Parametrize  the  straight  lines  thus  defined  in  the  f  =  1  plane  (these  are  just  the 
lines  along  the  vectors  (2.3b))  by  the  angle  9,  measured  positive  in  the  c»unterdockwise  sense 
from  the  x-axis,  with  the  straight  line  oriented  as  shown  in  Fig.  2.1  .  The  normal  (pointing  from 
the  b  side  to  the  a  side)  to  each  plane  is  given  by 

n  -  (n„  n„  «^.)  =    ±  (  S/,a,b)  tan  6  -  Sgia,b),  -tan  9,  1  )  ,  (2.4) 

where  the  +  sign  is  used  for  --^  s  9  ^  -^  and  the  -  sign  for  -^  ^  9  ^  -^  .   These  planes 
will  obey  the  entropy  condition  (1.3)  for  the  jump  a  -» Z>  provided 

£(*)  ~±(k-b)  hs^(a,b)  -  Sg(k,b)]  -  tan  9  [5/a,i)  -  S^k,b)]\  ^  0  ,  (2.5) 

for  every  k  such  that  fe  s  jfc  £a.    The  +  and  -  signs  are  correlated  with  the  angle  9  as  above. 
Thus 

Theorem  2.3:  A  point  of  jump  (t  =  l,x,y),  in  whose  neigbourhood  the  solution  is 
piecewise  smooth,  with  the  one-sided  limit  values  a  and  b  will  obey  the  entropy  condition  (1.3)  if 
and  only  if  its  tangent  vector  (2.3b)  is  oriented  so  as  to  obey  the  conditions  of  (2.5).  □ 

For  the  case/  ■  g  (2.5)  takes  on  a  simpler  form, 

E(k)  =  ±(k-  b)  [S^a,  b)  -  Sy(k,  *)]  (1  -  tan  9)  a  0  .  (2.6) 

The  next  theorem  follows  directly  from  (2.6). 

Theorem  2.4:  Let  (t  =  l,x,y)  be  a  point  of  jump  of  a  piecewise  smooth  solution  to  the 
Riemann  problem  for  the  case/  ■  g.  Then, 

a)  if  Sj(a,b)  &  S^k,b)  for  every  it  in  [6,  a]  the  point  of  jump  will  obey  the  entropy  condition 

(1.3)  if  and  only  if  its  tangent  vector  (2.3b)  (oriented  as  in  Rg.  2.1)  lies  in  the  angular  range 

^  s  9  s  -^  .    (Note  the  correlation  for  choice  of  '+'  and  '-'  sign  with  angle  as  given  by 

(2.4)  whidi  must  be  taken  into  account  to  obtain  this  result.) 

b)  if  S^a,b)  s  S^k,b)  for  every  k  in  [b,  a]  the  point  of  jump  will  obey  the  entropy  condition 
(1.3)  if  and  only  if  its  tangent  vector  (2.3b)  (oriented  as  in  Fig.  2.1)  lies  in  the  anguJar  range 


-^  s  e  ^  -J-  .   (Note  the  correlation  for  choice  of  '+'  and  '-'  sign  with  angle  as  given  by  (2.4) 

which  must  be  taken  into  aax)unt  to  obtain  this  result.) 

c)  if  the  states  a  and  b  are  such  that  the  quantity  S^a,b)  -  SXk,b)  changes  sign  for  some  k^ 
where  fc  <  *o  <  o ,  the  point  of  jump  cannot  obey  the  entropy  condition,  o 
Thus  for  /  ■  g ,  the  states  a  and  b  can  occur  as  the  left  and  right  limiting  values  of  u  at  a  point  of 
jump  that  obeys  the  entropy  condition  if  and  only  if  either  S^a,b)  ^  S^k,b)  or  SAa,b)  ^  SAk,b) 
for  every  *  in  [a,  b]. 

In  two  spatial  dimensions,  shock  waves  are  then  smooth  surfaces  of  discontinuity,  each  point 
of  which  is  a  point  of  jump  and  satisfies  the  entropy  condition  (1.3). 

The  RBi-efaction  Waves 

Let  (f,  X,  y)  be  a  reg\ilar  point  on  a  smooth  level  surface  of  a  piecewise  solution  u  to  a  two 
dimensional  Riemann  problem.  From  (1.1)  we  conclude  that 

[!,/'(«),  «'(«)] 
is  tangent  to  the  surface  at  t,  x,  y. 

However,  by  the  self-similarity  of  the  Riemann  solution,  the  vector 

[\,x/t,y/t] 
must  also  be  tangent  to  the  point.  Consequently,  the  vector 

[  0.  x/r -/'(«),  y/f -«'(«)]  (2.7) 

is  tangent  to  the  point  on  the  level  surface. 

Proposition  2.5:     a)  The  intersection  of  a  smooth  level  surface  with  the  plane  r  =  1  is  a 
straight  line  segment  defined  by  the  vector  (2.7). 
Let  u  =  a  and  u  =  b  denote  two  such  level  surfaces.  Then, 

b)  K  fl  ^  i,  the  two  level  surfaces  cannot  cross. 

c)  If  a  =  fc,  the  only  point  that  can  be  common  to  the  two  level  surfaces  and  the  r  =  1  plane  is 
X  =  fia),  y  =  g'{a). 
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Proof: 

a)  This  follows  immediately  from  (2.7).   The  level  surface  is  therefore  a  segment  of  a  plane  in 
t,  X,  y  space. 

b)  IS  a  ¥'  b,  then  the  value  of  u  on  the  curve  along  which  the  two  level  surfaces  cross  is  not 
unique,  contradicting  the  assumption  that  the  point  belongs  to  a  level  surface. 

c)  Follows  as  a  consequence  of  the  proof  in  part  b)  and  (2.7).  D 

A  rarefaction  wave  is  now  definable  in  a  manner  analogous  to  the  one  spatial  dimension  case 
(See  [1]).  Let  n  be  the  unit  normal  to  a  smooth  level  surface  labelled  u,  at  the  point 
p  =  (tp,  Xp,  yp).  By  the  smoothness  of  u,  there  exists  a  p  ±  8^  n  neighbourhood  of  level  surfaces 
for  sufficiently  small  Sp  >  0  .  For  some  -8^  <  c  <  8^,  let  m,  label  one  such  level  surface.  In 
the  plane  defined  by  the  vectors  n  and  the  time  axis  f  these  level  surfaces  appear  as  straight  lines 
segments.  If  the  lines  defined  by  these  segments  intersect  at  some  point  (r,,  J,,  y^  for  t^  <  t  for 
every  0  <  |c|  s  8^  then  the  level  surface  is  a  rarefaction  wave.  We  note  that  a  rarefaction  wave  is 
the  union  of  diaracteristic  lines. 

Proposition  2.5  states  that  two  rarefaction  waves  (of  the  same  value  «  =  r)  can  intersect  in 
the  f  =  1  plane  at  the  point  x  =  f'(r),y  =  g'{r).  This  is  a  reflection  of  Lemma  2.2  which 
characterizes  the  continuity  (but  possible  lack  of  differentiability)  of  the  solution  across  the  plane 
X  =  y.  This  situation  is  pictured  in  Fig.  2.2b  for  the  waves  of  two  rarefaction  fans  along  the  plane 
x  =  y. 

As  in  the  one  dimensional  case,  a  rarefaction  fan  is  defined  as  an  open  set,  all  points  of 
which  are  in  rarefaction  waves.  A  constant  state  is  a  domain  (connected  open  set)  in  space-time 
t,x,y  over  which  the  solution  is  constant. 

Explidt  Form  of  the  Non-Linear  Waves 

Based  on  the  results  of  §3  of  [1]  the  forms  of  the  non-linear  waves  for  the  case  f  ^  g  can  be 
explicitly  characterized.  The  form  of  a  typical  rarefaction  fan  is  illustrated  in  Rg.  2.2a  as  it 
^spears  in  the  plane  f  =  1. 

Fig.  2.2c  displays  the  3  general  forms  of  shock  waves.  The  rarefaction  fans  that  are  used  to 
distinguish  the  three  forms  are  drawn  in  as  dotted  lines.   The  shock  form  labelled  D^  separates 


-9- 

two  regions  of  constant  u  values,  denoted  u  =  a,  and  u  =  b.  The  form  labelled  X^^  separates  a 
region  of  constant  u  value  from  a  fan  of  rarefaction  lines.  The  F  shock  line  discussed  by  Wagner 
[3]  corresponds  to  a  special  case  of  this  type  of  shock  form.  The  shock  labelled  2^  separates  two 
fans  of  rarefaction  lines. 

The  shock  form  Z^  will  be  a  straight  line  segment.  The  line  defined  by  this  segment  will 
pass  through  the  point  (S^(a,b),  Sg(a,b))  in  the  jc,  y  plane. 

Fig.  2.3  shows  a  I^^  shock  separating  a  region  of  constant  u  =  a  from  a  rarefaction  fan 
terminating  at  its  upper  end  with  the  value  u  =  b.  For  clarity  in  the  figure,  we  represent  the 
rarefaction  fan  only  by  its  boundary  wave,  u  =  b  and  one  interior  wave.  We  first  obtain  the 
equation  of  the  shock  for  the  case  0  <  p(r)  <  it/2,  for  every  r  in  the  rarefaction  fan,  and  for 
it/2  <  a  <  tt.  We  parametrize  the  shock  by  the  u  values  of  the  rarefaction  fan.  Thus  the  shock 
plane  int,x,y  space  has  the  parametric  form 

r(a,ft)-(f,A:(r)f,y(r)r).  (2.8) 

A  normal  to  the  shock  surface  is 

«  =  (  ^'(r)y(r)  -  a:(r)/(r),  /(r),  -x'(r)  )  .  (2.9) 

Letting  >(r)  ■  ■y(r)  be  the  unknown  function  to  be  solved  for,  we  have  from  Fig.  2.3 

^(r)  =  fir)  +  (g'(r)  -  -y(r))  tan  P(r)  ,  (2.10) 

yir)  =  yir)  . 

The  jump  condition  (1.2)  applied  to  this  shock  gives  the  ordinary  differential  equation  (for  r  =^  a) 
for  -y(r), 

,.  .  _    ilir)  -  5,(r,fl))  [f"(r)  +  g"(r)  tan  P(r)  +  (g'(r)  -  y(r))  tan'P(r)] 
^^'^  /'(r)  -  5/r,a)  +  tan  P(r)  [g'C-)  -  Vr,«)]  '      ^        ^ 

having  boundary  condition 

Sf(a,b)  +  5,(fl,b)  tan  g  -  f'jb)  -  g'jb)  tan  p(r) 

yi")  - "^rr^ •  (2-12) 

'^  '^  tana  -  tan  3(r)  ^        ^ 

The  presence  of  the  non-linear  term  due  to  the  dependence  of  p  on  r  precludes  further  analysis 
without  information  concerning  the  form  of  p(r).    However,  if  we  assume  all  lines  in  the 
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rarefacdon  are  parallel,  v^cfa  in  practical  constructions  is  often  the  case,  ^  becomes  a  constant 
and  (2.11)  reduces  to 

,.  .  ^  (■y(r)-5,(r,a))(/"(r)-Kg"(r)tanP) 
^  ^'^        fir)  -  S^r,a)  +  (g'(r)  -  5,(r,a))  tan  P    '  ^  '     '' 

Equation  (2.13)  has  the  form 

-y'('-)  =  /(/•)  bir)  -  S^ir,a)]  , 

thus  defining  l(r).  The  solution  to  the  system  (2.12),  (2.13)  can  be  written 

r  1 

y'ir)  =  S^{r,a)  -h  (-/(*)  -  S^{b,a))  e"        -  /  |i-iA__Aj^J   ,   r         ^        (2.14) 

We  point  out  that  our  stated  restriction  on  the  angles  a  and  p  consfitutes  no  real  restriction. 
For  any  combination  of  angle  a  and  p  the  derivation  of  (2.11)  and  solution  of  (2.12),  (2.13) 
proceeds  in  an  analogous  manner.  For  some  angles  it  is  more  advantageous  to  let  x{r)  rather  than 
y{r)  be  the  unknown  function  to  be  solved  for.  The  cases  a,  p  =  n  tt/l  represent  special  cases  in 
which  the  equations  simplify. 

The  derivation  of  the  general  differential  equation  for  shocks  of  type  'L„  proceeds 
analogously  to  that  for  type  2^^  but  leads  to  a  much  more  unmanageable  equation.  We  therefore 
derive  the  equation  for  this  shock  type  for  the  case  in  which  the  lines  in  each  rarefaction  fan  are 
parallel  to  one  another  (see  Fig.  2.4). 

We  again  parametrize  the  shock  plane  by  the  «  values  of  one  of  the  two  rarefaction  fans  (say 
for  the  lower  fan  pictured  in  Fig.  2.4).  Equations  (2.8)  and  (2.9)  still  hold,  but  now 

_.,^  =   ie'(«(r))  -  g'(r)  -H  f'{u{r))  cot  a  -  /'(r)  cot  P 
'^^'  cot  a  -  cot  p 

^r^.  =  fiuir))  -  f'(r)  +  g'(u(r))  tan  a  -  f('(r)  tan  P  O  15) 

^^^  tan  a  -  tan  p  ^  '    ^ 

Substitution  of  (2.15)  in  (2.9)  and  application  of  the  jump  condition  (1.2)  gives  an  equation  for 
uir) 

u'(r)  [  g"(u(r))  tan  a  +  /"(«(r))  ]  [/'(r)  -  5/«(r),r)  +  (  g'(r)  -  5,(«(r),r)  )  tan  p  ] 
=  [  g"(r)  tan  p  +  /"(r)  ]  [/'("(O)  "  S^uir),r)  +  (  g'(u(r))  -  5,(«(r),r)  )  tan  ot  ]  , 
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(2.16) 
having  the  boundary  condition  u(a)  =  b  (see  Fig.  2.4)  which  will  be  a  known  value. 
For  the  case/  ■  g,  the  angles  drop  out  of  (2.16)  giving 

^  ^       /"("(r))         fir)  -  5>(r),r)  ^^^  ^ 

No  solution  is  possible  to  (2.16)  or  (2.17)  without  knowing  the  form  of  /.  If  we  assume, 
y(r)  =  /  r^,  the  solution  to  (2.17)  is  u(r)  =  -r  +  a  +  b  what  a  and  6  are  as  defined  in  Fig. 
2.4  .  Substitution  of  the  solution  to  (2.16)  or  (2.17)  in  (2.15)  gives  the  parametric  equation  for 
the  shock  curve. 

As  a  final  comment  we  note:  a  rarefaction  line  u  =  a  will  terminate  tangentially  on  a  shock 
(S^^orS^)  in  the  plane  r  =  1  at  a  point  p  where  the  limiting  u  values  on  each  side  of  the  shock 
wave  are  M  =  a  and  u  =  b  only  if  the  point  (f'{a),  g'(a))  coincides  with  the  point 
(5/a,fc),  S^(a,b)). 

In  the  half  spaces  x  <  y,  x  >  y,  irregular  points  in  the  solution  occur  as  intersection  curves 
of  shock  surfaces,  and  at  the  discontinuities  in  the  initial  data  in  the  plane  f  =  0.  In  the  f  =  1 
plane  these  irregular  points  are  therefore  points  at  which  N  (^  3)  shocks  meet.  The  line  x  =  y 
may  also  contain  irregular  points,  each  corresponding  to  a  rarefaction  wave  M  =  ainthex>>' 
half  plane  meeting  a  rarefaction  wave  u  =  a  inthex<>'  half  plane  continuously  but  not 
smoothly.  Irregular  points  corresponding  to  the  meeting  of  two  shock  waves  (one  from  each  half 
plane)  can  also  appear  on  this  line.  A  third  possibility  is  the  termination  of  a  shock  wave  on  the 
line  X  =  y.  For  this  to  occur  the  shock  wave  must  terminate  with  zero  strength. 

3.  Two  Dimension  Method  of  Solution  Construction  -  Examples 

Construction  of  solutions  to  2-dimensiona]  Riemann  problems  was  first  approached  by 
Guckenheimer  [2]  who  gave  the  solutions  to  two  example  problems.  Additionally,  he  listed 
several  of  the  principles  concerning  the  restriction  of  the  solution  in  the  plane  f  =  1.  Wagner  [3] 
has  constructed  solutions  to  (1.1)  for  the  case  that  the  functions  /  and  g  are  both  convex 
(analogously  both  concave)  and  that  the  initial  data  is  constant  in  the  four  quadrants  of  the  (x,  y) 
plane.   For  the  case  /  ■  g,  his  solutions  obey  the  entropy  condition  (1.3).   He  demonstrates  that 
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if  /  and  g  are  'sufficiently  close",  additional  restrictions  on  the  derivatives  of  the  functions/  and  g 
are  reqxiired  for  the  /  ■  ^  solutions  to  hold  for  the  /  #  g  case. 

Having  classified  the  most  general  non-linear  waves  that  can  appear,  the  construction  method 
these  two  papers  have  introduced  can  be  used  to  construct  the  entropy  obeying  solution  to  the 
general  2-d  Riemann  problem. 

The  method  of  solution  construction  is  the  foDowing:  The  construction  is  done  in  the  f  =  1 
plane.  The  curve  (/"'(«),  g'(u))  and  the  set  of  points  ((S^v.h'),  5^(v,w))  which  dictate  the 
orientation  of  rarefaction  lines  and  shock  curves  are  identified.  The  non-linear  waves  and  their 
analytic  forms,  and  the  irregular  points  from  which  the  solutions  will  be  composed  are  identified. 
The  angle  restrictions  which  the  entropy  condition  imposes  on  the  tangent  vectors  to  the  non- 
linear waves  are  obtained.  In  practice  this  imposes  limits  on  the  extent  to  which  a  particular  shock 
wave  form  may  appear  in  the  plane  in  a  particular  region  and  determines  the  'fitting  together'  of 
the  various  non-linear  elements.  The  solution  in  the  region  x^  -t-  y^  ^  a^  for  some  a  is  obtained 
as  the  solution  of  non-interacting  1-dimensional  Riemann  problems.  The  solution  is  then 
extrapolated  into  the  region  x^  +  y^  <  a^  by  'fitting  together'  non-linear  waves.  In  the  following 
sub-sections  we  display  the  entropy  obeying  solutions  for  a  few  classes  of  Riemann  problems  for 
the  cases/  ■  g. 

For  the  case  /  «  g  the  non-linear  wave  forms  and  the  irreguJar  points  of  the  solution  follow 
from  §2  and  [1].  This  analysis  also  holds  for  functions  /  and  g  which  differ  by  a  multiplicative 
constant  since  this  can  be  rotated  into  a  problon  /  ■  g.  We  note  that  example  1  of  [2]  falls  into 
this  category. 

3.1.  A  Convex  /  Exampk 

We  illustrate  the  case  for  /  having  no  inflection  points  with  initial  data  consisting  of  three 
wedges  of  constant  states  a ,  b  and  c  (a  <  b  <  c)  centered  on  the  origin  inthex^y  half  of  the 
plane.  Fig.  3.1  shows  the  six  unique  (ie.  'topologically  different')  entropy  obeying  solutions. 

By  Theorem  2.4(a)  it  follows  that  the  solutions  satisfy  the  entropy  condition.  By  Theorem 
2.4(a),  the  three  wedges  can  be  changed  in  size  relative  to  one  another  anywhere  in  the  range 

— J-  s  9  s  -J-  and  the  solutions  will  not  change  their  topological  form. 
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Completely  analogous  solutions  and  conclusions  foUow  from  Theorem  2.4(b)  for  the  concave 
case  of  /  ■  ^. 

Shocks  of  type  2^  do  not  appear  in  the  solutions  in  Fig.  3.1.  For  them  to  appear  in  the  case 
of  convex  /  it  is  necessary  to  go  to  initial  data  consisting  of  4  wedges  in  the  half  plane  x  ^  y. 
Twenty-four  'topologically  distinct'  solutions  exist,  one  of  which,  showing  a  shock  of  type  2^  we 
sketch  in  Fig.  3.2  . 

3.2.   A  Single  Inflection  Point  Example 

We  now  consider  the  solution  for  /  »  g  obeying  the  properties 

1)  /:  [a,  b]^[c,d]  ,     a  <  b,  c  <  d,     a,b,c,d   i    R  , 

2)  /  is  strictly  monotonic  on  [a,b]  , 

3)  /  has  a  single  inflection  point  at  «,  m(a,  b)  , 

4)  /"(«)  >  0  on  [a,  «,),  /"(u)  <  0  on  («,,  b]  . 

Without  loss  of  generality  we  shall  take  a  =  0,  6  =  1,  c  =  0,  d  =  1.  f  has  the  property  that 
for  every  u  €  [0,  1],  there  exists  a  unique  u'  d  [0,  1]  such  that /'(«")  -  Sy{u,  u").  The  meaning 
of  the  iterated  notation  u"  is  then  to  be  understood  as  («")'. 

Siocks  of  type  S^r  separating  a  fan  of  rarefaction  lines  of  values  u  >  u,  from  a  fan  of  values 
u  <  Uj  appear  in  the  solution  to  the  two  dimensional  Riemann  problem  for  this  /  function.  If  a 
rarefaction  wave  has  value  m  =  r  on  one  side  of  the  shock,  the  rarefaction  wave  leaving  on  the 
other  side  will  have  the  value  u  -  r'  (defined  above)  and  will  leave  the  shock  line  tangentially. 
The  form  derived  for  shocks  of  type  2^  in  §2  does  not  apply  since  there  it  was  assumed  that  the 
angles  a  and  P  of  the  two  rarefaction  fans  were  known;  consequently  the  shape  of  the  shock  in 
that  analysis  was  found  by  determining  the  state  values  of  the  rarefaction  lines  on  one  side  given 
the  state  values  on  the  other.  Here  knowledge  of  the  states  on  either  side  of  the  shock,  and  the 
implication  that  one  set  of  rarefactions  leave  tangentially  determines  the  shape  of  the  shock  curve. 
Therefore  we  determine  parametrically  the  shape  of  this  shock  type  and  show  the  shock  obtained 
obeys  the  entropy  condition  (1.3)  for  the  case/  «»  g. 
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Consider  the  Riemann  problem  shown  in  Rg.  3.3(a)  with  w  and  v  as  in  Fig.  3.3(b).  The 
solution  is  shown  in  Fig.  3.3(c)  where  the  shock  of  type  2^  is  labelled  ^(r,  r').  The  analysis  for 
the  form  of  "^  proceeds  as  in  §2.  Parametrize  ^  as  in  (2.8)  with 

x{r)  -  Ur)  ,  (3.1) 

y{r)  =  fir)  +  (f'(r)  -  ^r))  cot  a(r)  . 
implication  of  the  jump  condition  (1.2)  gives  the  equation  for  «{», 

.,(.  ^  f  f"(r)  a  +  cot  a(r))  +  (Mr)  -  /'(r))  cot'(«(r))  1  ,3  ^. 

(f'ir)  -  S^r,  r'))  (1  +  cot  a(r) 

with  boundary  condition 

^(,)  =   r(v)  (1  +  cot  gfv))  -  r(v')  fl  ^  cot  B(v'))  (3  3) 

cot  a(v)  -  cot  p(v") 

As  in  (2.11),  variation  of  a  along  the  rarefaction  introduces  a  non-linear  term  into  (3.2).  If  we 
assume  all  lines  in  this  rarefaction  are  parallel  (as  indeed  we  have  drawn  in  Rg.  3.3(c))  which  will 
be  the  case  in  the  solutions  we  display  in  this  p^jer,  (3.2)  simplifies  greatly  to 


(  »|;(r)  -  SAr,  r')  ) 
(/'(O-S/r.O) 


(3.4) 


having  the  solution 

r  r 

Mr)  =  S^r,  r')   +    [  ,|,(v)  -  5/v,  v")  ]  e'  -   J  e'  ^  'q,  ^  '        (3-5) 

where 

/(,)  =  f"{^)      ^    .  „  (3.6) 

/'(r)-S/f,0 

Note  that  given  r  in  [w*,  v],  the  form  of  the  function  /  guarantees  the  existence  of  r"  and 
for  every  k  such  that  r'  ■&  k  ■^  r  we  have  S^r,  r')  s  S^k,  r').  From  the  construction,  it  is  clear 
that  the  above  shock  obeys  the  conditions  of  Theorem  2.4(b)  and  hence  obeys  the  entropy 
condition. 


-15- 

We  return  briefly  to  Fig.  3.3(c)  to  discuss  the  other  curved  shock  which  is  composed  of  two 
shocks  of  type  2^,,  one  separating  the  fan  of  parallel  lines  m  =  w  -  «  =  v"  from  the  constant 
region  «  =  w,  the  other  separating  a  fan  of  non-parallel  lines  u  =  v'  -  u  =  w"  from  the  constant 
region  u  =  w.  The  two  shocks  have  the  same  tangent  line  where  they  meet.  The  construction 
and  Theorem  2.4(a)  guarantee  that  this  composite  shock  obeys  the  entropy  condition. 

In  Fig.  3.4  we  present  the  entropy  obeying  solutions  to  (1.1)  for  /  ■  g  with  /  as  described 
above  for  initial  data  of  three  wedges  of  constant  states  a,  b  and  c  {a  <  b  <  c)  centered  on  the 
origin.  We  consider  only  the  two  general  cases:  either  a<Mj<fc<c  or  a<ft<u,<c, 
otherwise  the  problem  reduces  to  the  purely  convex  or  concave  case  discussed  in  §3.1  .  By  the 
qualitative  symmetry  of  the  function  /  about  its  inflection  point,  there  are  only  six  qualitatively 
distinct  solutions.  As  previously  mentioned,  the  solutions  obtained  in  the  at  &  >>  half  of  the  plane 
can  be  obtained  independently  of  the  solutions  in  the  x  ^  y  half. 

4.  The  Two  Dimensional  Rlemann  Problem  for  Two  Phase  Flow  io  Porous  Media 

For  incompressible  two  phase  flow  in  a  porous  medium,  with  capillarity  effects  and  gravity 
ignored,^  the  system  of  equations  to  be  solved  in  a  source  free  region  is 

-^-Hv'(5,/')V/(5)  =  0,  (4.1a) 

Vv'(j,P)  =  0.  (4.1b) 

where  s  is  the  volume  fraction  (saturation)  of  one  of  the  two  phases  (the  other  fraction  being 
1  -  j)  and  P  is  the  pressure  field  in  the  jwrous  medium.  The  field  i?  is  the  total  fluid  (consisting 
of  both  phases)  velocity.  The  velocity  i7  is  generally  assumed  to  be  proportional  to  V/*  (Darcy's 
Law)  and  when  the  form  of  this  relation  is  specified  the  system  (4.1)  can  be  solved.  A  form 
commonly  used  ior  f{s)  is  (immiscible  flow  [6] ) 

s'^  +  r(l  -  sy 

where  r  is  the  ratio  of  the  two  phase  viscosities  (we  assume  constant  viscosities).   The  fractional 
We  have  also  set  the  rock  porosity  function  <{)(:?)  =  1  . 


-16- 

flow  curve  given  above  then  has  all  the  properties  of  the  /  function  described  in  §3.2  . 

One  method  of  solving  (4.1)  numerically  consists  of  solving  (4.1a)  and  (4.1b)  sequentially. 
Thus  for  purposes  of  solving  the  hyperbolic  equation  (4.1a)  the  velocity  field  can  be  taken  as  some 
given  vector  field  \'(i).  For  two  dimensional  problems,  the  hyperbolic  problem  is  therefore 

f  +  v,(x,  >)igl  +  v^,(x,  y)^  =  0  .  (4.3) 

Shocks  in  two  phase  systems  are  very  rapid  (in  this  model,  discontinuous)  dianges  in  the 
phase,  from  a  region  that  has  largely  phase  1  in  its  pore  spaces  to  that  having  largely  phase  2.  To 
be  specific  we  will  caU  the  phases  oil  and  water,  with  s  =  0  corresponding  to  pure  oil  and  5  =  1 
corresponding  to  pure  water.  The  system  (4.1)  describes  a  secondary  oil  recovery  procedure 
where  injected  water  is  used  to  force  reservoir  oil  towards  production  wells.  The  oil-water  shock 
surface  is  typically  called  a  'bank'  (water-bank  if  the  water  is  displacing  the  oil,  oU-bank  if  the  oil 
is  displacing  water).  We  can  therefore  study  the  dynamics  involved  in  the  interaction  of  two  such 
banks,  for  example  when  the  water  banks  from  two  separate  injection  wells  meet,  by 
{^Tproximating  the  problem  as  a  Riemann  problem. 

Fig  4.1(a)  depicts  the  problem  to  be  solved  and  Fig.  4.1(b)  its  approximation  as  a  Riemann 
problem.  Let  v  denote  a  value  of  v,  averaged  over  the  interaction  area;  similarly  let  w  be  an 
averaged  value  of  v  .  Equation  (4.3)  becomes 

ii  +  yM£i  +  ^M^  =  0 .  (4.4) 

dt  dx  dy 

This  is  now  of  the  form  (1.1)  with  /  and  g  differing  by  a  multiplicative  constant.  We  can 
therefore  solve  (4.4)  in  rotated  coordinates  J,  y  where  f  =  g. 

Fig.  4.1(c)  displays  the  solution  to  the  Riemann  problem  of  Fig.  4.1(b)  in  some  rotated  ^,  t] 
system  (the  rotation  of  course  depends  on  the  relative  strengths  of  v  and  w;  we  have  assumed 
V  >  0,  w  >  0).   Fig.  4.1(d)  displays  the  solution  when  a  and  c  are  interchanged. 

We  note  that  for  flow  with  gravity,  the  flux  function  associated  with  two  phase  flow  has  two 
inflection  points,  as  opposed  to  a  single  inflection  point  for  the  gravity  free  flow  case. 
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5.   CoDchiding  Remarks 

The  general  solution  to  the  2-dimensional  Riemann  problem  for  the  case  f\™ fi^ f, 
f  (.  C^  :R  -R,  f  having  at  most  one  inflection  point  has  been  shown  to  be  characterizable  in 
terms  of  non-linear  (rarefaction  and  shock)  waves  in  a  manner  analogous  to  the  1 -dimensional 
Riemann  problem.  We  believe  this  diaracterization  in  terms  of  non-linear  waves  is  the  correct 
formulation  in  which  to  approach  the  construction  of  solutions  to  2-dimensional  Riemann 
problems  for  the  general  case  /j  =^  /i.  We  conjecture  that  the  set  of  rarefaction  and  shock  waves 
described  in  this  paper  will  provide  the  complete  solutions  for  many  cases  /.  =^  Z,.  A  regularity 
theorem  giving  conditions  under  which  all  2-dimensional  Riemann  problem  solutions  are  piecewise 
smooth  would  immediately  give  greater  generality  of  appbcation  to  the  non-linear  wave  forms 
given  here. 

These  solutions  to  2-dimensional  Riemann  problems  also  supply  a  set  of  problems  for  testing 
of  finite  difference  schemes.  The  richness  of  structure  of  these  solutions  lends  itself  to  this 
purpose. 
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Flgure  Captions 

Fig.  2.1 

The  tangent  line  to  a  point  of  jump  on  a  shock  curve  in  the  r  =  1  plane.  At  the  point  the 
shock  separates  the  two  states  a  and  b,  b  <  a.  The  tangent  passes  through  the  point 
(5^a,  b),  Sgia,b))  at  angle  9  with  respect  to  the  positive  x-axis  with  the  orientation  that  the 
larger  state  a  is  always  to  the  left  of  the  shock  line  when  viewed  from  the  above  point 
towards  the  origin. 

Hg.  2.2 

a)  A  typical  rarefaction  fan  as  it  appears  in  the  f  =  1  plane,  b)  Two  rarefaction  fans 
meeting  continuously  but  not  smoothly  on  the  x  =  )»  line  in  the  f  =  1  plane,  c)  The  three 
general  forms  of  shock  waves  for  the  two  dimensional  Riemann  problem  under  assumptions 
discussed  in  the  tejrt.   Dotted  rarefaction  waves  are  drawn  to  clarify  the  picture. 

Fig.  2.3 

The  most  general  example  of  the  shock  type  I^^ .  For  purposes  of  clarity,  only  two  lines  of 
the  rarefaction  fan  that  the  shock  separates  from  the  constant  region  m  =  a  are  drawn  in. 
Fig.  2.4 

A  particular  example  of  shock  type  S^  where  the  two  rarefaction  fans  on  either  side  are 
composed  of  separately  parallel  lines.  For  darity,  only  two  of  the  lines  in  each  rarefaction 
have  been  drawn  in. 

Fig.  3.1 

The  six  qualitatively  unique  entropy  obeying  solutions  in  the  region  x  ^  y  to  (1.1)  for 
f\  ^  f\  ™  f,  f  having  no  inflection  points,  with  initial  data  composed  of  three  wedges  of 
constant  u  value  (a  <  b  <  c)  centered  on  (0,  0). 

Fig.  3.2 

A  sketch  of  a  four  wedge  initial  data  solution  to  (1.1)  showing  the  appearance  of  a  Z^ 
shock.  The  2^  shock  has  been  drawn  as  a  straight  line  which  would  be  the  case  for 
/(u)  =  u^  /  2.  Note  that  the  2^  shock  segment  joins  smoothly  to  the  2.^  segment  with  the 
same  slope.  The  value  a  can  be  determined  from  the  intersection  of  the  2^^  shock  and  the 
rarefaction  line  «  =  c.  In  this  example  the  states  are  ordered  d  <  c  <  b  <  a. 
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Fig.  3.3 

A  two  dimensional  Reimann  problem  for  the  function  /  described  in  §3.2.  (a)  The  initial 
data,  (b)  The  function/  defming  various  states  found  in  the  solution,  (c)  The  solution. 

Fig.  3.4 

The  siT  qualitatively  distinct  solutions  in  the  x  &  y  half  plane  for  the  three  wedge  Riemann 
problem  discussed  in  §3.2  Note  that  one  of  the  three  initial  data  states  a,  b  and  c  is 
separated  from  the  other  two  by  a,,  the  inflection  point  of/. 

Fig.  4.1 

(a)  The  interaction  of  two  water  banks  in  a  two  dimensional  oil  reservoir,  (b)  The  Riemann 
problem  approximation  to  the  interaction,  (c)  and  (d)  The  solutions  to  the  Riemann 
problem  approximation  for  two  example  situations. 
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